I would really love some input and feedback about what language we should be using; especially from people for whom English is not their first / primary language.

Moi! I'm one of those :) As a long-time programmer, I like dimension and dimensionality. This being said, I do know what rank means and it's certainly a good word if we're going to work with matrices. The problem that I see it that the complete definition of rank (at least according to google AI) is

The rank of a matrix is the maximum number of linearly independent rows or columns it has, representing the dimension of the vector space spanned by its rows (row space) or columns (column space).

(Emphasis mine) If this definition is right (and I think it is), then I don't think rank is a good name. A mathematician using ndarray might think that rank() returns/calculates the actual mathematical rank. IIUC, a 3x3 matrix of zeros has a rank equal to 1, no?

Yes, rank does have that meaning! I have been considering swapping our use of Rank to be the name of the generic type currently called NDim, and making NDim the name of the associated type of the trait. I'd also change the trait name from Ranked to Dimensioned, and change the function name from rank to ndim. Thoughts?

Yes, imo, it's a good idea.

Hi! I consider myself an applied mathematician (working in quantitative finance). I have a preference for using the word rank for denoting the number of axes of a multi-dimensional array; while not universal, this does seem common in numpy, R, MATLAB, and Julia communities.

Some other scattered thoughts:

• while the rank of a matrix has a well-defined meaning in mathematics, in the context of a multi-dimensional array library I think there is less chance of confusion. It might be better to call this concept/function matrix_rank() (as numpy does), as it applies only to 2-dimensional arrays
• dimensionality might be hard to remember. For a given array whose dimensions are (c1, c2, ..., cp), there are multiple notions that reasonably map to the colloquial notion of dimensionality:
  1. p (the number of dimensions)
  2. c1 * ... * cp (the dimension in Euclidean space in which the array may be embedded);
  3. c1 + ... + cp (if each coordinate had names, this would be total number of names captured by the array)
     whereas rank seems to have consensus to mean 32 bit vs. 64 bit indexing #1. Pr

As for the function name to return the "rank", my preference is ndim(). It agrees with numpy/R (ndim) and Julia (ndims).

Anyhow, just my two-cents. Cheers, and thanks for maintaining this library.